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ABSTRACT 

Recent theoretical claims not withstanding, muons can and do oscillate. Muons pro- 
duced in association with neutrinos (if masses and mixing angles are nonzero) exhibit a 
joint oscillating spatial distribution. The possible use of muon oscillations as a probe of 
neutrino mass and mixing parameters is discussed using very simple physical arguments. 
Space-time oscillations in the secondary decay vertices of muons (produced by pion decay) 
persist after summing over all undetected neutrinos. 



1. Introduction 

Classic precision measurements of {g — 2) have been made for electrons [1] and muons [2- 
5]. For the latter case, muons are injected into a ring with a uniform applied magnetic 
field B. The experimental quantum survival probability of the muon in the ring to decay, 
(via fx~ e~ +i'f>+T^n) with a detected electron energy above a threshold value, has been 
fit to the theoretical oscillating form 

P _ ^-Mc^Tt/E { ^ + A cos(Ot + <P) \ 

^M--e-+.,+.eW-e l + Acos(P J- ^ ' 

In Eq.(l.l), M, E, and F"^ represent (respectively) the mass, energy, and intrinsic lifetime 
of the muon; t is time in the laboratory reference frame, and O = (((7 — 2)eB /(2Mc)) is 
the experimental frequency measured in the muon survival probability due to quantum 
mechanical amplitude interference. Eq.(l.l) has been previously derived by us from the 
Dirac equation [6]. 

In view of the fact that the theoretical muon decay oscillations in Eq.(l.l) have been 
verified by some of the most precise experiments ever performed in high energy physics, it 
is sad to read a recent statement by Dolgov, Morozov, Okun and Schepkin[7] that "Muons 
do not oscillate." The seminal muon oscillation work by Russian chemists [8] at their own 
institute in Moscow was ignored. Some incorrect physics in the work of Dolgov, Morozov, 
Okun and Schepkin[7], including (but not limited to) the fact that their equations are not 
dimensionally correct, is summarized in the Appendix. 

Our purpose is to formulate the problem of neutrino induced muon oscillations using 
the notion of "quantum beats" which are discussed in recent text books[9-10] on quantum 
mechanics. The basic notion of a quantum mechanical decay exhibiting quantum beat 
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oscillations may be formulated as follows: (i) The survival amplitude at time t for a 
normalized state l^') is given by 

S{t) = {^\e-'"'/^\^). (1.2) 

(ii) The probability density for the state to have energy E is given by 

p{E) = {^\5{E-H)\^). (1.3) 

From Eqs.(1.2) and (1.3) it follows that the survival amplitude for a state is the Fourier 
transform of the energy probability density 

S{t) = J p{E)e-'^'/^dE, (1.4) 

yielding the survival probability 

P{t) = \S{t)\^ ^ J J p{E)p{E')e~'^^-^'^*/^dEdE' . (1.5) 

Eqs.(1.2)-(1.5) are exact. If the energy probability density p{E) may be written as a sum 
of peaked resonances, then it follows rigorously that the survival probability P{t) exhibits 
quantum beat oscillations in time. The proof follows directly from Eq.(1.5). 

Eq.(1.5) (in a fixed reference frame) and its more general Lorentz invariant forms of 
expression are the basis of the discussion which follows for the theory of neutrino mass 
mixing induced muon decay oscillations. In Sec. 2 quantum beat oscillations are discussed 
using the standard formalism successfully applied to atomic and molecular experimental 
decays. In Sec. 3, the Lorentz invariant version of the theory is developed. In Sec. 4 the 
four momentum distribution for the muon in the reaction vr"*" — > + is discussed. In 
a model for which the neutrino has three possible mass eigenstates, and in which there 
is mixing, the (incoherent) momentum distributions obtained by us are identical to those 
conventionally used by others to analyse previous experiments. In Sec. 5, the surival prob- 
ability of the muon is computed and the interference effects (which are present if the 
neutrinos have both mass and flavor mixing) are discussed. In the concluding Sec. 6, the 
rigorous relationship between the four momentum probability distribution and the survival 
probability is reviewed. An incoherent sum of the four momentum probability distribution 
may nevertheless exhibit quantum beat oscillations in the survival probability. 

2. Quantum Beat Oscillations in Decays 

Consider the case of quantum beats in a decay. Suppose that the initial state of the 
system under consideration has the form 

1^') =Ci|^'i) + C2|*2>, (2.1) 

with the assumption that l^'i) and |^'2) have no common decay modes. Explicitly, 

(^2|e-*^*/^|^'i> = ^(*2|n,owt>e-*-^"-°"**/'^(n,o'ut|*i> = 0, (2.2) 
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so that 

{^2\S{E-H)\^^) = 0, {^^\S{E-H)\^2) = 0. (2.3) 
Eqs.(1.3), (2.1) and (2.3) imply 

p{E) = |ci| V(^) + \c2fp2{E), (2.4) 

where 

Pi{E) = (*i|5(E - H)\^^), p2{E) = {^2\S{E - ff)|*2>. (2.5) 

Since the two states and {'^2) have no common decay modes, the energy prob- 
abihty distribution p{E) for the state 1^^) is a superposition of probabilities of the energy 
probability densities pi{E) of state l^^i) and P2{E) of state 1^^2)5 including the probabili- 
ties |c,f = |(*,|^)|'. 

The superposition of energy probability densities in Eq.(2.4) by no means implies that 
there is a superposition of survival probabilities in time. In fact, Eqs.(1.4) and (2.4) imply 
a superposition of survival amplitudes having the form 

S{t) = |ci|'5i(t)+ |c2|'52(t). (2.6) 
If the two states decay into different channels according to an exponential decay law 

5i(t) = e-riVSg-iSitA, S2{t) = e-r2*/2e-i£^2V?i, (2.7) 

then the survival probability 

Pit) = |5(t)|2 = |ci| V^i* + |c2| V^^^* + 2|ci|'|c2|'e-(r^i+r^^)*/2cos(a;i2t), (2.8) 
where 

nuJi2 = El - E2. (2.9) 

Eqs.(2.8) and (2.9) are fundamental equations for the quantum beat oscillations in 
decays which constitute (for example) commonplace experimental phenomena in atomic 
and molecular physics and chemistry. 

3. Lorentz Invariant Space-Time Oscillations 

The amplitude that a normalized state describing a physical system in the neigh- 
borhood of the space-time origin survives a displacement to the neighborhood of the space- 
time point X may be defined by 

S{x) = (^je^^-^/'^l*), (3.1) 

where 

P=(P,iy/c) (3.2) 

is the operator four momentum. Note that Eq.(3.1) reduces to Eq.(1.2) for decays in the 
center of mass frame[ll] for which the three momentum may be set to zero. 
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In a Lorentz covariant manner, one may start from Eq.(3.1) and write 

Six) = J a{k)e'''-''d% 



where the four momentum distribution is given by 



a{k) = 

If the initial state is in a superposition of finitely many states, 



(3.3) 



(3.4) 



(3.5) 



no two of which have common decay modes, then the four momenta probability densities 
are given by the (incoherent) superposition weighted by the probabilities |cjp 



where 



ajik) 



(3.6) 



(3.7) 



Eqs.(3.6) and (3.7) are the Lorentz invariant form of Eqs.(2.4) and (2.5). The survival 
amplitude follows from Eqs.(3.3) and (3.6) to be of the superposition of amplitudes form 



j 

Clearly, the survival probability exhibits quantum interference effects 

P{x) = \S{x)\'^ = J2 |cincj|^5*(a;)5j(a;). 



(3.8) 



(3.9) 



and quantum beat oscillations in decays have been described in a Lorentz invariant form. 

4. Muon Momentum Distributions for tt"*" —>//"'' + 

We consider a massive neutrino model with mixing and define the following parame- 
ters: For three flavors of neutrino (ve, i^^-, and v^.) there is a rotation matrix 



(4.1) 



where ni, n2, and denote neutrino mass states with masses mi, m2, and respectively. 
In the reaction 

7r+ ^ + zv^, (4.2) 
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the four momentum distributfon for the fx'^ has the form of Eq.(3.6); i.e. 
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where (Jmuon{k,mj) is the normalized four momentum distribution of the obtained if 
the original muon produced in the pion decay were to recoil against a neutrino of mass 

Note that the four momentum probability density for the muon in Eq.(4.3), is an inco- 
herent probability superposition. We have already summed over the unobserved neutrino 
final states, both fiavor and momentum. The incoherent four momentum distribution for 
the muon has long been understood (and employed) in the analysis of experimental pion 
decays as a probe of possible neutrino masses[12]. Unfortunately, the previous experimen- 
tal resolution for measured muon momentum distributions did not allow for the definitive 
detection of any finite neutrino masses. However, the resolution did allow for an improved 
determination[12] of the pion mass. 

5. Muon Survival Probability from 7v~^ ^ n'^ + 

From Eqs.(3.6), (3.8) and (4.3), the muon survival amplitude against e'^+P/j^+Ve 
from muons produced in a pion decay is given by 

3 ^ 

where 

^muoni^j'^j) — J '^muon{k,Tnj)e''^ k (5-2) 

is the survival amplitude of the //"'" obtained if the muon were to recoil against a neutrino 
of mass rrij. The survival amplitude Smuon{x,mj) is very well known; i.e. 

Smuonix^rrij) ^ e'^y^l^e-^V^I-^c^ (5 3) 

In the pion rest frame 'Pj = {pj , ^M^c^ + |pjP) , 



|p,.| _ ^{M^ - (^M + mjr) {Ml - (M^ - m,)2) 



2M, 



(5.4) 



In any reference frame, the survival probability for the muon follows from Eqs.(5.1) and 
(5.3) to be 

3 3 

Pmuon{x) ^ \Smuon{x)\ ~ e"^^''^/'' ^ ^ | | \R^i\ COs[{j)j - Pi) ■ X /%) . (5.5) 



The presence of quantum beat oscillations in the muon survival probability, see Eqs.(5.4) 
and (5.5), is the central result of this work. The incoherent sum in the four momen- 
tum probability distribution in Eq.(4.3) exhibits quantum beat oscillations in the survival 
probability in Eq.(5.5). 

6. Conclusion 

In the pion decay 7r+ — > + z/^, the muon later also decays via ^ e+ + z/^ + 
i^e- Assume that the neutrinos have both mass and flavor mixing. The four-momentum 
distribution CTmuonik) of the muon produced by the pion will be a probability superposition 
as in Eq.(4.3). This result is well known[12], and has already been employed in the analysis 
of pion decay experiments searching for neutrino masses. For the muon survival amplitude, 
the well known Eq.(4.3) yields Eq.(5.1) and for the survival probability Eq.(5.5). 

The interference effects are manifest in the survival probability Eq.(5.5). The sec- 
ondary vertex from — > e'^ + i>i^ + Ue, distributed via Pmuon{x), may provide a more 
sensitive probe of neutrino masses and mixing than does the four momentum distribution 
via CTmuonik)- Further details about oscillations may be found in our earlier works[13-16]. 



APPENDIX 

In his lectures on quantum field theory, T. D. Lee[17] states "An equation in physics, 
say A = B, must satisfy the requirement that the dimension of A be the same as that of 
B: [A] = [S] . This seemingly elementary requirement can serve the useful purpose of 
verifying the correctness of one's equations, especially after a long calculation." 

Almost all the central probability equations in the work of Dolgov, Morozov, Okun and 
Schepkin[7], hereafter abbreviated DM0S[7], are dimensionally inconsistent and thereby 
fail to satisfy T. D. Lee's elementary requirement for correct physical equations. 

In Eq.(9), DM0S[7] define a dimensionless joint muon neutrino probability 

P^^ (x^, X.) = Pi + + 2plPlcosicf>, - cl>2) (DMOSl) 

where (3ja are seen to be dimensionless in Eq.(ll) of DM0S[7]. In Eq.(25), DM0S[7] define 
a single muon probability 

^m^(^m) = E / Pt{^^.^^)d^^. {DM0S2). 

The quantity on the left hand side of Eq.(DM0S2) is alleged to be dimensionless (in fact 
unity) from Eq.(27) of DMOS. The right hand side of Eq.(DM0S2) has dimensions of 
\length^^ . By virtue of dimensional analysis alone, their result for the muon probability 
distribution is wrong. The integrations are simply incorrect. Errors evident from dimen- 
sional considerations are also present (for example) in Eq.(26) of DMOS [7]. 

The dimensional errors are symptomatic of a confusion in the starting wave function 
Eq.(6) of DM0S[7] 

^pAx^l^xAx^) = \l^> (e~*'^^cos^|i/i > ^e-^'l'^ sine\v2 >). {DM OSS) 
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Imagine a world in which one of the neutrinos has a large mass and the other neutrino 
has virtually zero mass. An energetic pion decays at CERN so that the muon "ket" which 
recoils from the heavy neutrino decays in Switzerland and the muon "ket" which recoils 
from the massless neutrino decays in France. In the DM0S[7] work the "ket" in France 
and the "ket" in Switzerland are identical, i.e. |/U > in Eq.(DM0S3). We prefer to think 
that < fJ.FrancelfJ' Switzerland >= 0. At any rate, without this physical notion there is more 
than a little difficulty in properly normalizing wave functions. 
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